Vortex transmutation 
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Using group theory arguments and numerical simulations, we demonstrate the possibility of chang- 
ing the vorticity or topological charge of an individual vortex by means of the action of a system 
possessing a discrete rotational symmetry of finite order. We establish on theoretical grounds a 
"transmutation pass rule" determining the conditions for this phenomenon to occur and numerically 
analize it in the context of two-dimensional optical lattices or, equivalently, in that of Bose-Einstein 
condensates in periodic potentials. 



Vortices are a physical phenomenon common to all 
complex waves. Defined by a phase singularity im- 
plying the vanishing of the wave amplitude their pres- 
ence is ubiquitous in physics where examples of vor- 
tices can be found in as diverse systems as quantized 
superfluids and superconductors, Bose-Einstein conden- 
sates (BEC's), nonlinear optical structures, or low dimen- 
sional condensed matter or particle systems (for a review 
see 0, 0)- The possibility of changing vortex proper- 
ties using periodic systems is a natural step based on 
the known example of the different behavior of electrons 
with or without the presence of a crystal. Like electrons, 
properties of vortices in a lattice have been shown to be 
qualitatively different than in a homogeneous medium. 
Vortices have been numerically predicted to exist in two- 
dimensional (2D) arrays of coupled waveguides in 
2D periodic dielectric media with Kerr nonlinearities — 
and, equivalently, in 2D BEC's with periodic potentials — 
0,0 and in photonic crystal fibers with defects They 
have been experimentally observed in optically-induced 
square photonic lattices [3, la] • In all cases the presence 
of the periodic medium has a strong influence on vortex 
features thus opening a door for their external manipu- 
lation. In this Letter we will show how manipulation by 
means of an external system owning discrete rotational 
symmetry can even affect the most intrinsic feature of a 
vortex, its vorticity or topological charge, leading to a 
phenomenon that we refer to as "vortex transmutation". 

Angular momentum is conserved in a nonlinear 
medium with 0(2) rotational symmetry in the x-y plane 
described by a first-order evolution equation of the type 
L(\(j)\)cf) = —id(f>/dz for the complex scalar field (f>. If 
we consider a solution with well-defined angular momen- 
tum I G Z (i.e., an eigenfunction of the angular momen- 
tum operator —id/ 89: tfii = e tW fi(r)) at a given axial 
point zq, evolution will preserve the value of I for all z. 
In a system possessing a discrete point-symmetry (de- 
scribed by the C n and C nv groups) angular momentum is 
no longer conserved. However, in this case one can define 
another quantity m G Z, the Bloch or pseudo-angular 
momentum, which is conserved during propagation 



The pseudo-angular momentum to plays then the role of 
I in a system with discrete rotational symmetry. From 
the group theory point of view, the angular and pseudo- 
angular momenta I and to are also the indices of the 2D 
irreducible representations of 0(2) and C„, respectively 
EHIlIl 113 Unlike I, the values of to are limited by the 
order of the point-symmetry group C n : |m| < n/2 |9tlll|. 

The appearance of this upper bound for the pseudo- 
angular momentum to opens the interesting question of 
determining the behavior of solutions propagating in a 
0(2) rotational invariant medium with well-defined an- 
gular momentum / after impinging a medium with dis- 
crete symmetry of finite order in which the value of I 
exceeds the upper bound for pseudo-angular momentum. 
This question can be analyzed in the light of group the- 
ory. Let us consider a wave propagating in a 0(2) non- 
linear medium corresponding to a solution (pi (not nec- 
essarily stationary) with well-defined angular momen- 
tum / launched into a second nonlinear medium char- 
acterized by the C n group. The surface separating the 
two media defines an 0(2)-C n interface that we locate 
at z = 0. We assume evolution is first order in z: 
L(\(j}\)(j) — —idqb/dz. Evolution in the second medium 
is thus fully determined by the initial condition 0/(0). 
The initial field 4>i{0) will excite a different representa- 
tion of the C n group depending on the value of I. Once 
this second wave is excited, it will propagate in the C n 
medium by preserving its representation — defined by its 
pseudo-angular momentum to. Let tp m be a function 
in the representation of C n characterized by the index 
to. Let us determine now what values of I are allowed 
by symmetry to produce a nonzero projection of <?!>z(0) 
onto (p m for a given value of to. The projection coeffi- 
cient is given by c ml = J R2 tp* m {r,6)(j)i(r,6;Q), Since <p m 
and (f>i belong to representations of C n and 0(2), respec- 
tively, they both properly transform under a discrete ro- 
tation of order n: <p m (r,9 + 2n/n) = e™( 2,r / n )^ m (r,0) 
and 4>i{r,6 + 27r/n) = e lZ ( 27r /™Vz0, 0). Thus, by per- 
forming the change of variable 8 — > + 2ir/n in the 
definition of c m i one arrives to the symmetry relation 
Cml = cxp[i(l — m)2ir / n]c m i . The c m i coefficient is then 
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Figure 1: Allowed values of the pseudo-angular momentum m 
in a d medium in terms of the angular momentum I of the 
field <f>i impinging the interface from an 0(2) medium. 

zero unless the following condition is fulfilled: 

I — m = kn (k £ Z), where \m\ < — . (1) 

The m representation of C n is thus excited by initial fields 
having angular momenta I = m,m ± n,m ± 2n, .... In 
a (3(2)-0(2) interface each representation of angular mo- 
mentum m in the second medium is excited by one, and 
only one, angular momentum component I arising from 
the first medium and verifying I — m. Contrarily, in 
a 0(2)-C n interface the symmetry restriction (JIJ implies 
that, due to the cutoff \m\ < n/2 in the C n medium, 
there are infinite angular momenta I that can excite a 
given representation of index m in the second medium. 
This can be clearly seen in Fig. 2] in which we represent 
the permitted values of m for different values of I for 
the particular case of a second medium with fourth-fold 
symmetry (n — 4) according to the angular momentum 
"pass rule" in Eq. . The effect of the cutoff (m < 2) is 
apparent in this representation. 

Thus, when the incident field carries an angular mo- 
mentum I that overcomes the limiting value for pseudo- 
angular momentum in the second medium it will ex- 
cite a wave that will propagate with different constant 
pseudo-angular momentum m given by the "pass rule" 
(Q. This result is valid for waves verifying an equation 
of the type L{\<f>\)<f> = —id(j)/dz, linear or nonlinear, sta- 
tionary or evolving. A particularly interesting situation 
is that in which the incident field is a vortex field of 
the 0(2) nonlinear medium. This vortex field <j>J is a 
stationary solution of the evolution equation with well- 
defined angular momentum I ^ 0: LQ(f>J\)(pJ = —pL(f>J- 
We consider here individual "canonical" vortices with 
a single phase singularity (i.e., with only one point in 
which <j)J = 0): 4>J(r,0,z) = e m f^(r)e-^ z 13\. The 
vorticity or topological charge of such solutions will be 
given by the circulation of its phase gradient around 
the singularity v = (l/2n) § V arg(^)dr, which equals 
angular momentum for canonical vortices v = I. On 
the other hand, the propagating wave 4> m with pseudo- 



angular momentum m excited by <j)J will evolve in the 
C n medium. There are different options for the asymp- 
totic states of (f> m when z — ► oo. One possibility is 
this wave asymptotically tends to an stationary solution 
(t>m Z — *° 4>m{ r -> Q> z ) = e " n6 9m(r, 6>)e _v z in the represen- 
tation of C n given by the conserved pseudo-angular mo- 
mentum m. If 0^ has a single phase singularity then it 
will have the structure of an individual canonical discrete- 
symmetry vortex, its vorticity or top ological charge v' 
being directly given by m: v' = m |ll| . The formation in 
the C n medium of an asymptotic stationary state in the 
form of a discrete-symmetry vortex is a dynamical issue 
that depends on the structural parameters of the second 
medium as well as in the characteristics of the input vor- 
tex field c/>J. If dynamics allows the stabilization of the 
discrete-symmetry vortex solution, the 0(2)-C n interface 
will realize the mapping of a 0(2) vortex with charge 
v = I (exceeding the limiting value for pseudo-angular 
momentum) into a C n vortex with charge v' — m 0, 
such that v' < v. The "pass rule" for pseudo- angular 
momentum (QJ becomes a "pass rule" relating input and 
output vorticities: 

v - v' = kn (k e Z), (2) 

where v' presents a cutoff in terms of n given by:|w'| < 
n/2 (even n) and \v'\ < (n - l)/2 (odd n) [HI. Note 
that m = n/2 solutions are not vortices but nodal or 
dipole-mode solitons ^J. We will refer to the process of 
mapping an individual vortex into another with different 
topological charge as "vortex transmutation". 

We will provide now a physical example of system in 
which the phenomenon of "vortex transmutation" takes 
place. It is an optical interface separating two 2D dielec- 
tric media with Kerr nonlinearity, these two media being 
a homogeneous medium and a 2D square optical lattice. 
This system is equivalent to a 2D BEC in which a peri- 
odic potential is abruptly switched on. They constitute 
an 0(2)-Ci interface given by the following equation: 

(V 2 t -V(x : y,z)+ 1 (z)\cf ) \ 2 )4=-i^, (3) 

az 

in which V t is the 2D gradient operator and V(x, y, z) = 
Vq + 9(z)(V~i(x,y) — Vq) where 8(z) is the step func- 
tion and V and Vi{x,y) = Vt (cos 2 (^a;) + cos 2 (^y)) 
(V\ is the potential strength and A is the lattice spa- 
tial period) define the refractive index profile of the 
homogeneous medium and of the 2D optical lattice: 
V Q = -(n 2 - 71q) and Vi(x) = -(n 2 (x) - n 2 ), n be- 
ing a reference refractive index introduced by the slowly- 
varying envelope approximation. The nonlinear function 
7(z) = 7 + (1 — j)6(z) permits the nonlinear response 
of the system to be different in the two media. All dis- 
tances appearing in Eq.© are normalized and dimen- 
sionless (x = fcox', z = k^z'). In order to solve the 
evolution problem in this system we solve first Eq. J3J for 
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z < 0, which becomes an ordinary Nonlinear Schrodinger 
equation (NLSE) for a homogeneous medium. Since our 
aim is to evidence the phenomenon of "vortex transmuta- 
tion" we are interested in finding canonical vortex solitons 
of different charges in the homogeneous 0(2) medium: 
</^(x, z) — e lW fl(r)e~' lfJ - z . This can be done by standard 
methods. At a given value of I, a family of 0(2) vortices 
are found characterized by their power Pi = J R2 \4>i\ 2 and 
their propagation constant fj,, which are related through 
the relation Pi{n). In the case of a Kerr nonlinearity, fi 
behaves as a scaling parameter and Pi is /x-independent 
P|. Once the vortex solution (f)J is found, it is taken as 
an initial solution for propagating it in the 2D optical 
lattice (z > 0): 0(x,O) = <^(x,0) = e m fj{r). Thus we 
solve Eq.© for z > 0, which becomes a NLSE which the 
periodic potential V\ (x) , with the previous initial condi- 
tion. This is solved numerically using a standard split- 
step Fourier evolution method. 

According to our previous symmetry arguments, the 
evolution of the </> wave for z > has to occur in a way 
that the "pass rule" for angular momentum JQ) is ful- 
filled. The 0(2) vortex soliton <j)J carrying angular mo- 
mentum I will excite a propagating wave 4>m for z > in 
a representation of C4 with pseudo-angular momentum m 
given by Fig. ^ Indeed, numerical evidence of this "pass 
rule" is obtained by analyzing the rotational symmetry 
of the evolving field. By construction, the input momen- 
tum is / since we choose the solution to be of the form 
(j)J(x, z) = e iW ff{r )e"^ z for z < 0. In order to check the 
symmetry properties of the solution for z > 0, we numeri- 
cally evaluate the rotated field <j)(r, 9, z) = <fi(r, 9 + tt/2, z) 
at every step in z and compare it to its unrotated value 
(j)(r,9,z). If 4> belongs to the m representation of C4, 
4>{r, 9 + tt/2, z) = e im7r/2 0(r, 9, z) and the ratio (j>/<f> will 
have a constant value for all x 6 M 2 (with the excep- 
tion of x = 0, where rotations are ill-defined) and z > 0: 
<t>/<j) = e im7r / 2 . If this condition is satisfied the value of m 
can be directly extracted from the numerical ratio (f>/4>. 
Indeed, the independence of the 4>/<j) ratio from trans- 
verse coordinates is numerically verified at every axial 
step, which permits to evaluate m for different values of 
z > 0. Results are shown in Fig.0 These results nicely 
confirm the general condition JQl and, more specifically, 
they satisfy the graphical rule represented in Fig. ^ f° r 
an 0(2)-C/t interface. 

Once the angular momentum "pass rule" is checked 
there persist the question of the fake of the propagat- 
ing wave in the C4 medium. As predicted by theory, m 
is numerically conserved during evolution. However, the 
asymptotic behavior of the 4> m evolving field can be very 
different depending on the parameters of the incident 
vortex field (its power P and its propagation constant 
/j,) and of the characteristics of the periodic potential 
Vi(x) (the potential strength V\ and the lattice period 
A). Our interest lie in obtaining asymptotic stationary 
states which can be described as individual or canonical 
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Figure 2: Numerical confirmation of the angular momentum 
"pass rule" for different values of the angular momentum I of 
the incident vortex field: (a) I = 1 and m = 1; (b) I = 3 and 
m = —1; (c) I = 5 and m = 1; (d) I = 7 and m = —1. 








Figure 3: Amplitudes and phases of input 0(2) vortices and 
output C4 vortices for different values of the input vorticity v: 
(a) v = 3 and v = -1 (P = 2.5, fi = 2, Vi = 2, A = 1.3); (b) 
v = 5 and v' = 1 (P = 3.0, fi = 3, Vi = 3, A = 2.5); (c) v = 7 
and v' = -l(P = 3.5, fj, = 3, Vi = 3, A = 3.5). 



discrete-symmetry vortices. This condition implies that 
the asymptotic field has to present a single phase singu- 
larity. In other words, we want to exclude multi- vortex 
or cluster excitations. In order to achieve this feature, 
we enlarge the optical lattice (by increasing its period A) 
according to the size of the input vortex for increasing 
values of I. Thus, in our simulations A is fixed by I. 

After performing many different simulations, we have 
indeed found numerical evidence of the "vortex trans- 
mutation" phenomenon. By playing with the input pa- 
rameters P and fi and the lattice strength V\ and pe- 
riod A, we have been able to find asymptotic station- 
ary states 0^ = e lm0 g m {r 1 9)e~ 1 ^ z for different values 
of the input vorticity value v = I. The vorticity of the 
output field can be only v' = ±1 because of the vortic- 
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ity cutoff for a C4 system (recall that m = ±2 solutions 
are not vortices but nodal or dipole-mode solitons 
In Fig. |31 we show the amplitudes and phases of input 
and output vortices for different input vorticity values 
v. All of them verify the vorticity "pass rule" ©■ The 
"vortex transmutation" phenomenon only occurs when 
|w| > 2. Similar results are found for the corresponding 
input anti- vortices with negative values of v. When, for 
fixed v = I (fixed A), the election of P, fi, and V\ is not 
adequate the asymptotic solution can be non-stationary. 
We observe two different scenarios besides the station- 
ary regime: discrete diffraction of the input wave in the 
optical lattice and self-focussing instability leading to fil- 
amentation of the field. A thorough analysis of multiple 
configurations permits to elaborate a "vortex transmu- 
tation" phase diagram where the three different regimes 
can be recognized. As an example, in the phase dia- 
gram shown in Fig. 0] we observe the "vortex transmuta- 
tion" region (shaded) differentiated from the diffraction 
(white) and self-focussing instability (light shaded) re- 
gions as a function of the power P and propagation con- 
stant /i of the input vortex field at fixed V\ . Analogous 
phase diagrams are found for different values of V\ . It is 
interesting to analyze the evolution of different "vortex- 
transmuting" configurations by monitoring the evolution 
of the power P'(z) and average propagation constant 
fi'(z) = J <fi*(-id/dz)4>/ J 4>*4> (defined both on the fi- 
nite domain of the numerical solution) in the C4 medium. 
These quantities are z-dependent, in general. However, 
they become independent of z when we analyze a station- 
ary solution; thus we expect (P'(z), fi'(z)) (P', //) 
for asymptotic stationary states. Every input 0(2) vor- 
tex characterized by the initial values (P, /j.) defines then 
a different trajectory in the P'-// plane. In Fig. 0] (inset) 
we show four different trajectories mapping 0(2) v = 3 
vortices with different (P, ^1) initial values into asymp- 
totic C4 vortices with charge v' = — 1 characterized by 
their (P', //) values. It can be checked numerically that 
these values lie on the same P'(p') curve found in Ref.|l4| 
for stationary vortices with charge w' = -1 in an iden- 
tical square optical lattice. By launching a whole family 
of initial vortices we have been able to asymptotically 
reproduce the entire P'(p') curve of C4 vortices. It is re- 
markable that the asymptotic C4 vortices in the Fig. 0] 
inset have been checked to be stable under small per- 
turbations whereas the original 0(2) ones are not 
0. Thus the "vortex transmutation" phenomenon not 
only permits to change the charge of unstable input vor- 
tices but it can also help to transform them into stable 
structures. Inversion of the vortex charge has been ob- 
served in "noncanonical" vortices in free-space and occurs 
through its dynamic propagation |l3l |. Here, however, all 
vortices involved are "canonical", the key point to the 
"vortex transmutation" phenomenon to occur being the 
suitable matching between angular and pseudo-angular 
momentum at the 0(2)-C n interface. Despite our model 
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Figure 4: "Vortex transmutation" phase diagram for v = 3 
vortices in an optical lattice with Vi = 2. Inset: Evolution in 
the {P 1 , n') plane of four initial v = 3 vortices characterized by 
different initial (P, /i) values. The solid line curve corresponds 
to the stationary v' = — 1 vortices as in Ref . ["Ml] . 



refers to an specific system, the theory of transmuting 
vortices is general and applies to any system given by an 
equation of the type L{\<f)\)<f) — —id<j)/dz in the presence 
of an 0(2)-C n interface. Hence, it is expected this phe- 
nomenon to occur in a wide variety of physical systems 
as those mentioned in the introduction of this Letter. 
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